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Abstract 

Frequent lane-changes in highway merging, diverging, and weaving areas could disrupt traffic 
flow and, even worse, lead to accidents. In this paper, we propose a simple model for studying 
bottleneck effects of lane-changing traffic and aggregate traffic dynamics of a roadway with 
lane-changing areas. Based on the observation that, when changing its lane, a vehicle affects 
traffic on both its current and target lanes, we propose to capture such lateral interactions by 
introducing a new lane-changing intensity variable. With a modified fundamental diagram, we 
are able to study the impacts of lane-changing traffic on overall traffic flow. In addition, the 
corresponding traffic dynamics can be described with a simple kinematic wave model. For a 
location-dependent lane-changing intensity variable, we discuss kinematic wave solutions of the 
Riemann problem of the new model and introduce a supply-demand method for its numerical 
solutions. With both theoretical and empirical analysis, we demonstrate that lane-changes could 
have significant bottleneck effects on overall traffic flow. In the future, we will be interested in 
studying lane-changing intensities for different road geometries, locations, on-ramp/off-ramp 

flows, as well as traffic conditions. The new modeling framework could be helpful for developing 
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ramp metering and other lane management strategies to mitigate the bottleneck effects of lane- 
changes. 

Keywords: LWR model, fundamental diagram, lane-changing intensity variable, Riemann 
problem, kinematic waves, supply-demand method, bottleneck effects, NGSIM data 


1 Introduction 


An understanding of the evolution of traffic dynamics is the foundation of 


analysis, managemen t, control, and planning. In the celebrated LWR model l|Lighthill and Whitham 


1955 


Richards. 


transportation network 


1956f) . vehicular traffic is viewed as a continuous fluid flow, and traffic dynamics are 


described by the changes in {x, t)-space of three aggregate variables: density p, speed v, and flow-rate 
q. The model is can be written as 

dp dpV{p) 


dt 


dx 


= 0 , 


( 1 ) 


which is based on three assumptions: 


• The fundamental law of traffic flow: 


q = pv, 


( 2 ) 


Traffic conservation 


^ + ^ = 0 

dt dx 


(3) 


Fundamental diagram of speed-density relationship 


V = V{p). 


(4) 


The first two assumptions can be derived from the continuum assumption of traffic flow (IHaberman 


1977L Chapters 59, 60). An equilibrium speed-density relationship assumption (|1|) has been shown 
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to be valid in stationary states ( 

Del Castillo and Benitez 

1995 

). Sue 

h a relationship can also be 

derived from various car-following models in steady states ( 

Gazis et ah. 

1961 

) and ( 

Haberman 

1977 


Chapter 61). In this sense, the speed-density relationship in dH) captures the longitudinal interactions 
among vehicles, and the trad itional LWR model is consistent with car-following behaviors at the 


aggregate level (INewell 


196111 . In the LWR mod el, traffic dynamics are described by combinations 


1974 1. and the LWR model is thus called the kinematic 


of shock and rarefaction waves (|Whithaml . 
wave model of traffic flow. 

Major roadways, however, usnally have multiple lanes, and vehicles can not only change 
speeds, but also change lanes. That is, vehicles can have both longitudinal and lateral move¬ 
ments. A lane-changing area is a region, where one or more streams of vehicles systematically 
change their lanes. Such areas can be near a merging junction and lane-drops, upstream to a di¬ 


verging junction, inside a weaving section^ or around a clover le af interchange ( Milam and Choa 


1998; 


Cassidy and Rudianakanoknad 


and accidents (e.g. 


Golob et al 


200511 . Since bottlenecks (IHall and Agvemang-Duahl 


199111 


200411 tend to occur in these areas, it is important to understand 


phenomena associated with lane-changing traffic. In this study, we assume that different lanes are 
balanced on average; i.e., vehicle speeds are the same across different lanes at the same location 
and time. In reality, it has been observed that traffic is nearly balanced inside major weaves, where 
speed difference fo r weaving and non-weaving vehicles is not statistically significant (about 5 mph) 


(Roessetal 


197411 . Actually, imbalance among different lanes is usually a reason for lane-changing. 


and lane-changing traffic can have balancing effect; i.e., under certain situations, lane-changes could 
smooth out differences between lanes, and the balancing effect could be beneficial to the whole traffic 
system in achieving higher capacity. 

At the micro scopic level, vehic les’ speed adjustment behaviors have been modeled by car¬ 


following models (jCazis et al 


196111 : similarly, lane-choice behaviors on why, when, where, and 


3 


















































how a vehicle changes its lane ar e mod eled by lane-changing models (|Gippi 


Toledo et al 


2003 


Kesting et al. 


1986 


Yang . 


1997 


20071). These models can describe detailed lane-changing be¬ 


haviors, but usually contain a large number of parameters and cannot provide intuitive descrip¬ 
tions of system-level effects of lane-changing traffic. At the macroscopic level, many studies have 
been carried out to understand various characteristics of lane-changing traffic, including exchange of 


flows between lanes (|Michalopoulos et al. 


1984; 


200311 . density oscillation and instability issues (IGazis et al 


Holland and Woods 


degree of First-In-First-Out violation among vehicles (IJin et al 


1962 


1997 




2002 


Munial and Pipes 


1971 


Goifman 


and the 


2006ll . In particular, characteris¬ 


tics of weaving sections have been extensively studied since the publication of Highway Gapacity 
Manual in 1950, and many methods have been proposed to analyze levels of service at we aving 
areas, usually measured by speeds of weaving and non-weaving vehicles (e.g. 


Leisch 


1979). In 


these methods, parameters include weave type, geometric configurations, the number of lanes, the 


(jGassidv et al 


1989; 


Gassidv and May 


1991; 

Ostrom et al.. 

1993; 

Windover and Mav 


199J), the 


distributi on of weaving and non-weaving traffic on rightmost lanes was also studied for different 


locations. 


Eads et al 


20oi 


1 200011 proposed a framework for evaluating dynamic capacities of a weaving 


section. In (jLaval and Daganzo 


20061 1. a hybrid model of lane-changing traffic was proposed, and 


it was shown that systematical lane-changes could cause a capacity-drop consistent with observa¬ 
tions. All these studies, however, do not provide a simple approach for analyzing the impacts of 
lane-changing traffic and the corresponding traffic dynamics at the aggregate level. 

In this paper, we attempt to fill this gap and develop a simple kinematic wave model of lane¬ 
changing traffic. Based on the observation that, when changing its lane, a vehicle is the leader of 
two following vehicles on its current and target lanes, we modify the speed-density relationship in 
dH) by adding an effective additional density, equal to total density of all lanes times a lane-changing 
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intensity variable (or simply intensity) e{x,t). With the new fundamental diagram for both car¬ 
following and lane-changing traffic, we then introduce a kinematic wave model based on the LWR 
model. When the lane-changing intensity is location-dependent, we analyze the new kinematic wave 
model as a system of hyperbolic conservation laws. We also calibrate a relationship between lane¬ 
changing intensity and traffic density with observed vehicle trajectories in a weaving section. Note 
that lane-changing intensity e{x, t) is determined by drivers’ lane-changing choices and characteristics 
in a road section during a time interval. Therefore, e{x, t) is highly related to the location in a lane¬ 
changing section, on-ramp and off-ramp flows, and roads’ geometric configurations. In this paper, 
we demonstrate that, once e{x, t) is known, we are able to evaluate the impacts of lane-changing 
traffic at the aggregate level within the framework of kinematic waves. 

The rest of the paper is organized as follows. In Section 2, we introduce a lane-changing 
intensity variable, a modified fundamental diagram, and a kinematic wave model for lane-changing 
traffic. In Section 3, we analyze the kinematic wave model as a nonlinear resonant system and 
propose a corresponding supply-demand method for numerical solutions. In Section 4, we calibrate 
a relationship between the lane-changing intensity variable and traffic density for a weaving section 
with NGSIM data. Finally, some implications and extensions of this study are discussed. 

2 A kinematic wave model of lane-changing traffic 

2.1 Model derivation 

Assume the longitudinal axis is a;, and the lateral axis is y. Then the movement of a lane-changing 
vehicle can be represented by a trajectory in (x, y, t) space. As shown in Figure [U during a time 
period of At, a lane-changing vehicle can move Ax and Ay in the x and y directions, respectively. We 
can see that the lateral displacement threshold of the lane-changing vehicle, or simply lane-changing 
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Figure 1: An illustration of a lane-change 

threshold, Ay, should be at least the width of the vehicle w. To accomplish a lane-change, a vehicle 
first usually signals its intention. In very sparse traffic, it can easily find a gap big enough to switch 
to its target lane without waiting too long or interrupting traffic stream on either its current or 
target lanes. However, in congested traffic, it usually has to slow down first and wait for a gap 
or speed up to squeeze in. In this process, following vehicles on its current lane have to maneuver 
accordingly, and following vehicles on the target lane may slow down or switch to other lanes to 
let it in. That is, a lane-changing vehicle has longitudinal influence on its current lane and lateral 
influence on its subjec t lane. In a lane-changing area, individual vehicles’ lane-changing choices 


may vary significantly (|GiDDsl . 


198d l. but our experience and many studies suggest that systematic 


lane-changi ng vehicles may disr upt traffic flow. 


Since (jGreenshields . 


193511 . speed- or flow-rate-density relations, v = V(p) and q = pV{p) 


respectively, have been used to capture drivers’ response to traffic environment at the aggregate 
level. Note that, traffic density p is traditionally defined as the number of vehicles per unit length of 
road. Thus traditional fundamental diagram captures only longitudinal interactions between vehicles 
when they follow each other, but not lateral interactions when they change lanes. Therefore, the 
traditional LWR model cannot be directly applied to describe traffic dynamics in lane-changing 
areas. 


Since a vehicle uses both its current and target lanes during a lane-change, we propose to double 
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the contribution of a lane-changing vehicle to total density. I.e., lane-changing traffic causes effective 
additional density. If denoting a lane-changing intensity variable, the effective total traffic 

density is given by 

p{x,t) = p{x,t){l + e{x,t)). (5) 

Furthermore, for lane-changing traffic, we use the following modified speed-density relationship: 

V = V{p) = V{p{l + e)). (6) 

Then, the fundamental diagram with lane-changing effect is 

q = pV{{l + e)p). (7) 

Correspondingly, the fundamental diagram without lane-changing effect is 

q = pV{p) = {l + e)pV{{l + €)p). (8) 

By comparing Q and we will be able to understand the impact of lane-changing traffic. 

With the fundamental diagram incorporating lane-changing intensity, we can model lane¬ 
changing traffic dynamics in the framework of kinematic wave theories as follows: 

Pt + {pV{p{l + e)))^ = 0, (9) 

which can be considered as an extension of the LWR model for equilibrium lane-changing traffic. 


2.2 Determination of lane-changing intensity in uniform traffic 


Generally e{x,t) is time- and location-dependent, and e{x,t) > 0 due to the bottleneck effect. Since 
^ = V'p^ w hich is negative in general, this modified speed-density relationship is consistent with 


observations (|Fazio and Rouphail 


1986f) that the total number of lane-shifts in a weaving section. 


proportional to e, is negatively correlated to weaving and non-weaving speed. At the microscopic 
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level, lane-changing intensity e(x, t) is determined by drivers’ lane-changing choices and character¬ 
istics in a road section during a time interval. At the macroscopic level, the lane-changing intensity 
variable in ([5]) are related to locations, road geometric configurations, on-ramp and ofF-ramp flows, 
and other traffic conditions. Therefore, in practical applications, e(x, t) has to be calibrated for 
different locations of a lane-changing area and traffic conditions. 

Here we attempt to determine lane-changing intensity when traffic conditions are uniform in 
the lane-changing region; i.e., traffic density is the same across the region and for all lanes, and all 
vehicles travel at the same speed. This is not meant to provide a generic formula for computing 
e in ®. Instead, we intend to show which parameters could affect the lane-changing effect at the 
aggregate level. 

For a lane-changing section of length L, we have the following quantities: p^c is the density 
of lane-changing traffic, pnlc the density of non-lane-changing traffic, p = Plc + Pnlc the total 
density, v the speed of both lane-changing and non-lane-changing vehicles, qlc = Plc'u the flow-rate 
of lane-changing traffic, and qNLC = Pnlcv the flow-rate of non-lane-changing traffic. Then, the 
time for all vehicles to traverse the lane-changing region is T = L/v. In order to understand how 
much effect lane-changing traffic can have on total traffic, we have to know the time for finishing a 
lane-changing maneuver and the number of lane-changes. 

We denote the total number of lane-changes in the lane-changing area during a period of T 


ditions and road geometric 

configurations ( 

Worrall and Bullen. 

197(1: 

Worrall et ah. 

197(1: 

Pahl 

1972; 

Chang and Kao, 

1991; 

Klar and Wegener, 

1999 

). Here, we are mostly interested in the re- 


lationship between lane-changing traffic flow and the number of lane-changes for relatively con¬ 


gested traffic. In ([Fazio and Rouphail . 


1986 ). it has been calibrated that the total number of 


lane-changes are linear to weaving, i.e., lane-changing, flows for different lane configurations of 











































weaving areas. In (|Fitzpatrick and Nowlinl . Il996r) . it is observed that a weaving vehicle generally 
takes 1.33 lane-changes in one-sided weaving operations on one-way frontage roads. Thus here we 
simply assume a linear relationship between the number of lane-changes and lane-changing flow; 
i.e., Nlc = oiqLcT = apLcL, where the coefficient a, the average number of lane-changes of each 
lane-changing vehicle, could be related to lane-configurations, number of lanes, traffic conditions, 
drivers’ ten dency to lane-changing, and the length of the sectio n. In addition, a could also depend 


on location (|Cassidv and May 


1991 


Windover and May 


1994 b 


We denote lane-change duration by tiCi which starts when a vehicle signals its lane-changing 
intention and ends when it finishes a lane-change. If the width of a lane is H, and the average 
lane-changing angle 9, we then have 


tbC = 


D 


!tan9 


( 10 ) 


Note that 0 is related to drivers’ behavior, road geometric configurations, and tra ffic conditions. Usu¬ 


ally t he time to complete a lane-changing man euver is around 2.5 sec on freeways (jWang and Prevedouroi 


19981 1. while about 10 sec on surface streets (jSheu and Ritchie 


200lll . 


Therefore, by doubling the number of lane-changing vehicles during their lane-changing periods, 
the effective number of vehicles at any moment is 


- pLT -b NLcthC j j tbC 

N = - - - = pL + apLcL—- 


Then we can obtain the lane-changing intensity variable, e, in ([5]) as 


NhCiLC NLctLC PLCtLC 
e = -— = a- 


( 11 ) 


pLT 


NT 


PT 


( 12 ) 


In this sense, the lane-changing intensity can be considered as ratio of the total lane-changing time 
to the total traveling time during a time interval on a road section. As expected, lane-changing 
intensity variable is determined by the number of lane-changes, the lane-changing duration, and 


9 






















traffic density in the lane-changing area. In particular, a and p^c or can depend on road 


geometry, locations, and traffic conditions. 


speed v=60 mph 



Figure 2: A lane-changing area downstream to a merging junction 


As a simple example, we consider a uniform, three-lane lane-changing region, which is down¬ 
stream to a merging junction connecting two three-lane roads as shown in Figure [2][jJ. The length 
of the lane-changing area is L = 1000 ft, the width of a lane D = 12 ft, and both branches have 
the same density of p/2; i.e., traffic density on each lane is p/6 before merging. We assume that all 
lanes in the lane-changing region are fully balanced; i.e., traffic density on each lane in this region 
is twice as that before merging, p/3. Therefore, one third of vehicles from the merging branch, i.e., 
p/6, have to stay on the rightmost lane in the lane-changing region with no lane-change, one third 
will change one lane to the middle lane, and one third will change two lanes to the leftmost lane. In 
this case, the total density in the lane-changing region is p, and the density of lane-changing traffic 
is that from the merging branch; i.e., pLC = p/3. Then we have the average number of lane-changes 
of a merging vehicle 

p/6 • 0 -b p/6 • I -b p/6 • 2 


p/3 


= 1.5. 


If we assume that vehicle speed v =60 mph and the lane-change duration sec, then the 

^Here we simplify the problem by assuming that all lane-changes occur in the region downstream to the merging 
junction. In reality, some vehicles might change their lanes upstream to the merging junction. 
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time for a vehicle to traverse the lane-changing region is T = Ljv=\\A sec, and we can find from 
m that the lane-changing angle 9 = 3.1 degrees. Then from (fT^ . we can have 


NLctLC = + ^L2tLC = 

6 6 2 


and the lane-changing intensity variable in (HU is 


1 2.5 

2TT4 


0 . 11 . 


Here we assume that lane-changes are evenly distributed over the whole road section. In reality, 
however, the distribution may be uneven, and we may have a location-dependent lane-changing 
intensity variable e{x). 


2.3 Bottleneck effects of lane-changing traffic 


We consider the following triangular fundamental diagram (jMunial et al.l Il97lt iHabermanl [1977 


Newell 


199311 . 


V{p) = 


Q{p) = 




0 < P< Pc', 


Pc . Pj—P 


Pc<p<Pj, 


vfP, 


0 < p < Pc', 


-p), pc<p<pj. 


where Vf is the free flow speed, pj the jam density, and pc the critical density where flow-rate, q = pv, 
attains its maximum, i.e. the capacity. The values of these parameters ar e: Vf = 65 mph, p, = 240 


vpmp l. Pc = Pj/d) = 40 vpmpl, and the capacity is Qmax = 2600 vphpl ([Del Castillo and Benitez 


19951) . After introducing lane-changing effect, we then have the following fundamental diagram. 


V{e,p) = 


Vf, 


0<p<pc/(l + e); 


(13) 


Pc/(l + e) <p<pj/(l-be). 


Pi 
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Q(e,p) 


0 < p < pc/(l + e); 


(14) 


Vfp, 

Pc/(l + e) <P<Pj/(l + e). 




Figure 3: The triangular fundamental diagram with constant lane-changing effect 


For one example, we assume constant lane-changing effect, e = 0.1. The relationships be¬ 
tween the standardized density, speed, and flow-rate are shown in Figured From the figure, we 
can see that: (i) lateral interactions can be omitted when traffic is relative sparse; (ii) there ex- 
i sts capacity reduction of 1 — 1/1.1=9.1%, which is consistent in magnitude with observations in 


(jCassidv and Bertini 


1999 


and (iii) the observed jam density is lower than maximum jam den - 


sity, which could be one reason of observing different jam densities (IDel Castillo and Benitez 


199511 . 


That is, for the same road section, constant lane-changes can yield a lower throughput. This clearly 
demonstrates the bottleneck effect of lane-changing traffic. 

For another example, we consider the following density-dependent lane-changing effect, which 


^Note that reasons of bottleneck were not discussed in this study, and we suspect that systematic lane-changes 
caused by merging traffic from on-ramps could have significant contributions. 
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has a jump at critical density, 


e(p) = 


0, p< Pc] 

Pi 


(15) 


15+2 p/pj ’ Pc- 

As shown in Figure|31 the lane-changing effect is appro ximately linear for congested traffic, and the re¬ 


sulted fundamental diagram has a shape of reverse-A (|Koshi et al 


19831 1 ■ Since lane-changing inten¬ 


sity is related to the location, road geometry, on-ramp and off-ramp flows, and other exogenous con¬ 


dition s, it is possible to have a discontinuous p — e relationship. In (ICassidv and Rudianakanoknad 


20051 1. for example, it was observed that lane-changing intensity is significantly higher when a 16- 


vehicle queue appears on the shoulder lane. Therefore, this discontinuous p — e relationship could 
also be a cause of capacity-drop: when traffic density increases over a critical value, the number of 
lane-changes significantly increases, and the overall capacity drops as a result. 



p/Pi 



Figure 4: A fundamental diagram of reverse-A shape with density-dependent lane-changing effect 
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3 Traffic dynamics for location-dependent lane-changing in¬ 


tensities 


In this section, we consider a special case, where the lane-changing intensity variable, e(x,t), is 
independent of traffic conditions and only a function of location denoted by e{x). Although simplified, 
this assumption does capture the most pronounced feature that lane-changing intensity is location- 
dependent. For this scenario, we can analyze ([5]) and understand some fundamental traffic dynamics 
in lane-changing areas. In the literature, there have been many methods for analyzing such an 


inhomogeneous LWR model with discontinuous flu x functions (refer to (Jin et al 


In this study, we follow an approach proposed in (I Jin and Zhang 


200311 . 


2 OO 9 II for a survey). 


3.1 A system of hyperbolic conservation laws 

For a location dependent lane-changing intensity e(x), we obtain a simple conservation law of lane¬ 
changing coefficient, 

et = 0, (16) 


which is equivalent to saying that e is time-indepen dent. 


Without considering inhomogeneities of links (|Jin and Zhane , 


2003f) or merging and diverging 


effect, we can have the following system of conservation laws, 


Ut+F{U), = 0 , 


(17) 


where U = (e, p), F{U) = (0, Q(e, p)), x G R, and t > 0. Note that the domains of e and p are 
e G [0,oo) and p G [0, respectively. 

Here we use a differentiable, concave fundamental diagram in our analysis. For example, we 
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can use the following maximum sensitivity model ( Del Castillo and Benite: 


dtel 


1995) 


V{p) 


Vf <1 — exp 


l-exp( ^(^-1) 

Vf p 


(18) 


where Cj is the shock wave speed for jammed traffic. In this fundamental diagram, V' < Q and 

< 0 , 

For the kinematic wave model of lane-changing traffic, dnD, the second term of flux, F{U)x, 
can be linearized as dF{U)Ux with 


dF 


0 0 

p^V V + p(lFt)V' 


whose two eigenvalues, or wave speeds, are Xo{U) = 0, and Ai(t/) = V + peV' 
corresponding eigenvectors are 


Rq — 

V + p{l + e)V' 

, i?i — 

0 


-p^V 


1 


Q'((l + e)p)l!l The 


Since V >Q and V < 0, it is possible that Ai = 0 = Aq. Therefore, (El) is a non-strictly hyperbolic 
system. 

In addition, if defining the critical traffic state t/* = (e*,p*) by Ai(t/*) = 0, we then have the 
following results. First, 


^Ai(D*) = =(l + e)Q"((l + e)p)|c,. <0, (19) 

since Q is strictly concave in (1 -|- e)p. Second, 


= pVV. < 0. 

From (HU), the wave speed Ai is decreasing with respect to traffic density. In addition, Ai(e,p = 
0) = Vf > 0 and Xi{e, p) = epV < 0 when (1 -|- e)p = pj. Therefore, for any weaving factor e*, 

®By assuming p = p/{l + e), we have that Q = Y^pV{p) = j^Q(p) and Xi{U) = Therefore, 

Ai(F) = g = Q'((l-t6)p). 
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we can find a unique critical state C/* = (e*,/9*). Further, since d\i/dp < 0, Q{U^) is indeed the 


maximum flow-rate for e* and equals ma,x{Q{e, p)\e = 0}, where max{(5(e, p)|e = 0} is the 
capacity without lane-changes, and reflects capacity-drop caused by lane-changes. Here we call 
p* = r(e,) as the transitional curve, since traffic states on its left side are under-critical (UC) while 
those on its right side are over-critical (OC). From these properties, (ITTl) is a nonlinear resonant 
system and, in the ne i ghbor hood of a critical state, the Riemann problem can be uniquely solved 


(jlsaacson and Temple 


1992|). 


3.2 Solutions of the Riemann problem 

We consider the Riemann problem of (El) with the following initial conditions, 

U{x,t = 0) = 


Ul, 

X < 0, 

Ur, 

X > 0, 


where Ul = [^l^Pl) and Ur = {^RiPr)- Since, in the Godunov method (| Godunov 


19591) . general 


initial conditions can be approximated by piece-wisely linear functions, we can solve the Riemann 
problem at each bo undary to find the flux and t hen update traffic conditions for the next time step. 


According to (llsaacson and Temole 


1992ll . the Riemann problem is solved by a combination 


of two basic types of waves, which are associated with the two eigenvalues of dF'. we have 0- or 
standing waves with wave speed Aq = 0, and 1-waves with wave speed Ai. Note that 0-waves 
carry contact discontinuities, and 1-waves include traditional shock and rarefaction waves. In the 
p — e phase plane, {Ul, Ur) yielding 0-waves forms a 0-curve, and similarly we can define a 1-curve. 
Further, 0-curves are the integral curves of eigenvector Rq, and 1-curves are those of Ri. They are 
pV{p{l + e)) =const and e =const, respectively. A 0-curve, a 1-curve, and the transi tional curv e are 


shown in Figure [5l For a Riemann problem, i ts solu tions must satisfy both Lax’s (|Lax . 


Isaacson and Temple’s (jlsaacson and Temple 


19721) and 


1992h entropy conditions, such that resulted waves 
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Figure 5: A 0-curve, 1-curve, and the transitional curve in (p, e)-plane 


increase their wave speeds from left to right, and 0-curves do not cross the transitional curve. Since 
both 0- and 1-waves are linear in the sense that U{x, t) = U{x/t), then from the Riemann solutions, 
we can obtain the boundary flux through x = Q as q = Q{U{x = > Q)). 

As shown in Figure El when is UC, i.e. to the left of the transitional curve in the phase 
plane, we have pL < r(eL) and can divide the p — e phase plane into four regions by the thick 
curves and obtain the four types of wave solutions for different values of Ur, which satisfy the two 
aforementioned entropy conditions. Here kinematic wave solutions of Type 1 are carefully explained, 
and solutions of Types 2 to 4 are given in Appendix. 

Type 1 When Ur is in Region 1, pR < T{eR), Q{Ur) < Q{Ul), and Q{eR,T{eR)) > Q{Ul). That 
is, Ur is UC, Q{Ur) is not greater than Q{Ul), but the capacity at cr is not smaller than 
Q(Ul). The Riemann problem is solved by a combination of a standing wave and a forward 
traveling rarefaction wave, with an intermediate state Ui as shown by the four figures of Figure 
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Figure 6: Four types of Riemann solutions when Ul is UC 

[3 The top left figure in the p — e phase plane shows that {Ul,Ui) is a standing wave, and 
the two states share the same flow-rate Q{Ul)- The bottom left figure shows two fundamental 
diagrams for when e = and e = respectively and the Riemann solution of a standing 
wave and a rarefaction wave. The top right figure shows characteristic waves in the x — t plane, 
and we can see the discontinuity caused by the standing wave at x = 0 and the rarefaction 
wave on its right side. The bottom right figure shows the profile of density at a time instant 
tf). In this case, since the flow-rate remains the same across x = 0, then the boundary flux 

q = Q{Ul)- 

As shown in Figure [51 when is OC, i.e. to the right of the transitional curve in the phase 
plane, we have pL > r(eL) and can divide the p — e phase plane into six regions by the thick 
curves and obtain the six types of wave solutions for different values oi Ur, which satisfy the two 
aforementioned entropy conditions. The solutions of Types 5 to 10 are also given in Appendix. 
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Figure 7: An example for wave solutions to Riemann problem of Type 1 


1996 


If we introduce new defini tions of local traffic supply and demand (jPaganzo 


Nelson and Kumar 


1995 


Lebacaue 


2004 1 as follows, 


S{e, p) = max{Q(e, p) : p < p < Pj}, 
D{e,p) = max{(5(e, p) : 0 < p < p}. 


( 20 ) 


This is equivalent to saying that 


and 


S{e,p) 


Q(e,r(e)), pe[0,r(e)], 
Q(e,p), pe(r(e),pj]; 


pe[o,r(e)], 

D{e,p) = I 

[ Q(e,r(e)), pe(r(e),Pj], 

where (5(e,r(e)) is the capacity at e. These definitions lead to that both local supply and demand 
are not smaller than flow-rate; i.e., S{e,p) > Q{e,p), and D[t, p) > Q{e,p). 
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Type 

Ul 

D{Ul) 

Ur 

SiUR) 

q 

1 

UC 

Q{Ul) 

UC, Q{Ur) < Q{Ul) < Q{eR,T{tR)) 

g(eK>r(efl)) 

Q{Ul) 

2 



Q{Ur) > Q{Ul) 

> Q{Ur) 

Q{Ul) 

3 



OC, Q{Ur) < Q{Ul) 

Q{Ur) 

Q{Ur) 

4 



UC, Q{tR,T{eR)) < Q{Ul) 

Q{f-R, r(e_R)) 

Qi^^R, r(e_R)) 

5 

oc 

Q(eL,r(eL)) 

OC, Q{Ur) < Q{Ul) 

Q{Ur) 

Q{Ur) 

6 


> Q{Ul) 

OC, Q{Ul) < Q{Ur) < Q(eL,r(ei)) 

Q{Ur) 

Q{Ur) 

7 



UC, Q{Ur) < g(eL,r(eL)) < g(eK,r(efl)) 

Qi^R, r(e_R)) 

g(eL,r(eL)) 

8 



Q{Ur) > g(eL,r(ez,)) 

> Q{Ur) 

g(eL, r(eL)) 

9 



UC, Q{tR,T{eR)) <Q{Ul) 

Qi^R, r(e_R)) 

Q{f-R, r(e_R)) 

10 



UC, Q{Ul) < Q{eR,T{eR)) < g(eL,r(ez,)) 

g(eK>r(efl)) 

Q{f-R, r(e_R)) 


Table 1: Boundary flux, downstream supply, and upstream demand for ten types of Riemann prob¬ 
lems 


20 

























































Figure 8: Six types of Riemann solutions when Ul is OC 

In Table [U we summarize the solutions of boundary flux for ten types of Riemann problems 
discussed in the preceding subsection and compute the corresponding upstream demand and down¬ 
stream supply functions. From the table, we can conclude that the boundary flux can be computed 
by the minimum of upstream demand and downstream supply; i.e., 

q = min{S'(C/i{),D(C/L)}. (21) 

Although (121(1 is proved for time-independent and differential fundamental diagram, we expect that 
the supply-demand dehnitions in (1201) can also be applied to the triangular fundamental diagram and 
to the scenarios when e is time-varying. Therefore, this method can be considered as a simplification 
and extension to the Riemann solver in the preceding subsection and is more efficient for numerical 
simulations. 
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4 Calibration of lane-changing intensity in a weaving section 


In this section, we consider lane-changing traffic on a multi-lane freeway section on interstate 80 
in Emeryville (San Francisco), California, as shown in Figure ID The freeway section has six lanes, 
where lane 1 is a car-pool lane, and an on-ramp from Powell Street and an off-ramp to Ashby Ave. 
The width of lanes 1 to 6 is 12 ft. Traffic direction is from south to north. The road section is covered 
by seven cameras, numbered 1 to 7 from so uth to north, and vehicle trajectories w ere transcribed 


from video data by FHWA’s NGSIM project ( FedgnlHighyg y ^^ministration 


are four data sets (ICambridge Svstematics. Inc.l . 


2004, 


2005ab 


2006 1. In total there 


Data set 1 contains trajectories 


of all vehicles every fifteenth second on December 3, 2003 between 2:35pm and 3:05pm; Data sets 
2 to 4 contain the locations of each vehicle every tenth second on April 13, 2005 between 4pm and 
4:15pm, between 5pm and 5:15pm, and between 5:15pm and 5:30pm, respectively. Note that the 
lengths of the study location in the first and the other three data sets are different: The Ashby Ave 
off-ramp is included in the first data set, but not in the other three. 


camera 1 camera 2 camera 3 camera 4 eamera 5 camera 6 

550 ft 470 ft 210 ft 170 ft 180 ft 270 ft 


camera 7 
1100 ft 


-300 ft-4 4- 


\--data ^'et 1: 2950 ft- 
data sets 2x3,4: 1750 ft--'- ^ ->1 

-^ — ’<1 - ^- 







Figure 9: 1-80 study area 

In Figure ITOl we demonstrate the trajectories in both x-y and t-y spaces of vehicle 576 in data 
set 2, when it switches from lane 6 to lane 5. Here y = 60 is the lane separation line, and the width 
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Figure 10: Trajectories of a lane-changing vehicle: (a) 'm x — y space; (b) int — y space 

of the vehicle is w =6.3 ft. In both figures, we demonstrate Acc and At for two different Ay, w or 
1.5rc. From both figures, we can see that the lane-changing angles is decreased by more than half 
and the lane-changing time is more than doubled when we increase Ay from w to 1.5?c. Note that 
TSic is still smaller than the width of a lane, 12 ft. From the figures, we can see that, after entering 
the target lane, the vehicle first stay near the lane separation line and then adjust to the center. 
Such behavior can be confirmed with video data by camera 6. 

Following the example in Section 2.2, we can theoretically derive the lane-changing intensity 
variable in (jl2ll and corresponding fundamental diagrams in ([8]) and ([7]) for a lane-changing area with 
a width oi L = Xb — Xa and time period with a duration of T = h —ta, shown in Figure [9l Assume 
that all lane-changes are induced by on-ramp vehicles and finished in a region of L downstream 
to the merge and that traffic flow is uniform with speed v and flow-rate q, respectively. Following 
the example in Section 12.21 we can estimate the number of lane-changes in L during T = L/v as 
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2.?)qonT, since one sixth of the on-ramp vehicles will make zero up to five lane-changes. Thus the 
lane-changing intensity variable can be estimated as 


e 


NhCtLC _ ‘2.5qonTtLC _ o c:„ 
pLT ~ 


( 22 ) 


Since qon is usually controlled by ramp meters and tiC is almost constant, e is generally not a 
constant in p or v. it is decreasing in p and increasing in v, when traffic is congested. In particular, 
if qon = 800 vph, p = 200 vpm, L = 900 ft, and t^c = 5 s, we have e = 8.15% for u = 60 mph. Note 
that (1^^ is just an estimation for the weaving region in Figure^ since (1) lane-changes can also be 
caused by vehicles leaving to the off-ramp or vehicles from or to other adjacent on- and off-ramps, 
(2) lane-changes induced by the on-ramp vehicles can occur upstream to the merge, and (3) fewer 
than one-sixth of the on-ramp vehicles can make five lane-changes, since the first lane is a car-pool 
lane. 


In the following, we consider a road section of 900 ft downstream to the merge in Figure [9] with 
Xa = 950 ft and Xb =1850 ft. The length is L =900 ft or 274 m. The time interval is [ti,ti+T], where 
the smallest ti is the time when the first vehicles pass xt, and the greatest ti + T is the time when 
the last vehicles pass Xa- For four data sets, r =39, 63, 86 , 100 s, respectively. Here we consider 


different lane-changing thresholds Ay. Refer to (j.lin 


20101 1 for detailed methods for computing 9, e. 


p, q, and V. 

With a lane-changing threshold of Ay = w, we have the following results. From Figure [TTl we 
find a linear relationship between density and lane-changing angle 


9 = -0.5016+ 0.0121p, 


with R-square 0.9441. Here we can also calibrate a relationship similar to (l22l) as 


e 


-0.0247 


24.2712 

P 
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with R-square 0.8077. Hereafter we only calibrate the exponential function, since it usually yields 
better R-square. From Figure [TH we can see that q = pV{p) has a capacity of 13408 mph when 
p = 233 vpm, and q = pV{{l e)p) has a capacity of 12369 mph when p = 215 vpm. Lane-changes 
can cause 7.75% reduction in capacity. 


(a) 



Figure 11: Characteristics of lane-changes at different densities with a threshold of w 

With a lane-changing threshold of Ay = 1.5?c, from Figure [131 find a linear relationship 
between density and lane-changing angle 

e = -0.3222 -h 0.0086p, 

with R-square 0.9139. We can also find the following relationship between density and e 

e = 0.5579e-°°°^®P, 
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Figure 12: Fundamental diagrams with a threshold of w: Blue dots with lane-changing effect, and 
red dots without lane-changing effect 


(b) 


1 

1 

+ % 4+ 



1 





% ° 

d> 0 » + 

^-1- 

0 + 

0+ = ^ 

OH 

0 + 

o-b + “^ 



_1 

_1 


_ 

_ 

_1 

_ 

_ 

_ 


o + 
o + 

-^0 


n. 


with R-square 0.8917. From Figure [141 we can see that q = pV{p) has a capacity of 15014 mph 
when p = 251 vpm, and q = pV{{l + e)p) has a capacity of 12369 mph when p = 215 vpm. The 
lane-changes cause 17.62% capacity reduction. Potential capacity drop is more than twice than that 
for Ay = w. This is consistent with the observation in Figure 1101 since t^c is more than doubled 
with Ay = l.Sw. 

For the whole weaving location, we have density-depende nt lane-ch anging intensity variables. 


More detailed analysis show that e is also location-dependent (Ijin 


2010r) . Note that the kinematic 


wave model ([5]) is developed for a roadway with both lane-changing and normal sections. If we can 
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(a) 



Figure 13: Characteristics of lane-changes at different densities with a threshold of 


obtain an e — p relation everywhere, m can be re-written as 

Pt + {pV{p{l + e{x,p))))^ = 0. (23) 

If we define Q{x, p) = pV (p(l -I- e{x, p))), this model is equivalent to an inhomogeneous LWR model. 
In this case, if we are able to introduce an inhomogeneity factor and write the equation above as a 
system of hyperbolic conservation laws, the analysis in Section 3 will apply for this general model. 
When we are not able to introduce an inhomogeneity factor, we can use the supply-demand method 
in mu for numerical solutions of the kinematic wave model of lane-changing traffic ®. 
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Figure 14: Fundamental diagrams with a threshold of |w: Blue dots with lane-changing effect, and 
red dots without lane-changing effect 

5 Conclusion 

In this paper, we presented a model for lane-changing traffic dynamics in the framework of kine¬ 
matic waves. A key assumption is that the disruption lane-changing effect can be modeled by a 
modified speed-density relation with a new lane-changing intensity variable. Then, for location- 
dependent lane-changing intensities, we studies kinematic waves arising in lane-changing traffic and 
corresponding new dehnitions of local traffic supply and demand for computing the flow passing 
through a boundary. With both theoretical derivations and observed data, we demonstrated that 
lane-changing intensities are highly related to road geometry, location, on-ramp/off-ramp traffic. 
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lane-changing time, vehicle speeds, and other traffic conditions. In particular, with 75 minutes 
of vehicle trajectories collected for a freeway section on interstate 80, we calibrated lane-changing 
intensities and corresponding fundamental diagrams for different lane-changing thresholds. The re¬ 
sults support a functional relationship between density and the lane-changing effect for the whole 
weaving section, and the lane-changing angles are highly related to traffic density. Furthermore, it 
was suggested that lane-changes could cause 8% to 18% reduction in capacity, depending on the 
definition of a lane-changing threshold. 

By incorporating lane-changing effects in the fundamental diagram, this study provides a simple 
framework to look into lane-changing traffic dynamics in the framework of kinematic wave theory. 
From this study, we can see that lane-changing traffic could cause capacity drop, different observed 
jam densities, and fundamental diagram of reverse-A shape. Lane-changing traffic can also affect 
the formation and dissipation of shock and rarefaction waves on a roadway. This simple model can 
be easily integrated into a commodity-based kinem atic wave simulation model of network traffic in 


(j.Tin and Zhang 


2004 


Jin and .Tavakrishnan . 


2005^ to study system-wide traffic dynamics. 


The model studied here bears certain limitations by assuming full balance among lanes and 
constant e. Thus our model is not intended for ramp weaves (HCM 2000), but for usually balanced 
areas near a lane-drop or merging junction. By assuming an equilibrium speed-density relationship, 
do]) is not able to described detailed traffic dynamics when vehicles are accelerating or decelerating. 
In this study, we only studied the bottleneck effects of lane-changing traffic. In relatively sparse 
traffic, however, lane-changes could actually benefit the overall traffic flow due to their balancing 
effects, and such effects have yet to be included in e. In addition, the new model considers the 
lateral interactions between vehicles in lane-changing areas at the aggregate level, and more detailed 
analysis of microscopic lane-changing decisions and maneuvers could give us more insights on the 
aggregate effects of lane-changing traffic. 
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With observed data, we calibrated a simple relationship between lane-changing intensity and 
traffic density. As shown in lane-changing intensity is also related to on-ramp and off-ramp 
flows. For the weaving section in Figure[9l however, we do not have off-ramp flows in data sets 2 to 4 
and therefore cannot study such relationships. In addition, it was shown that Ay could significantly 
impact lane-changing intensities and therefore capacity reductions, and it can be calibrated by 
comparing road capacities with lane-changes and those without lane-changes. For example, for the 
weaving section studied in Section 4, capacity reductions are about 8% and 18% when Ay is set to be 
a vehicles width and 1.5 times of a vehicles width, respectively. In the future, we will be interested in 
calibrating the lane-changing threshold with more data sets or more detailed analysis. In the future, 
we will also be interested in studying lane-changing intensities for different lane-changing areas, 
on-ramp/off-ramp flows, and other factors. Such a study would be helpful for understanding traffic 
macroscopic traffic dynamics in lane-changing areas and for developing possible ramp-metering and 
lane management strategies for improving the overall traffic flow. 
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Appendix: Kinematic wave solutions of Types 2 to 10 

Type 2 When Ur is in Region 2, Q{Ur) > Q{Ul). The Riemann problem is solved by a combination 
of a standing wave and a forward traveling shock wave, with an intermediate state t/i as shown 
in Figure HSl In this case, q = Q{Ul)- 



Figure 15: An example for wave solutions to Riemann problem of Type 2 
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Type 3 When Ur is in Region 3, pR > r(e_R) and Q{Ur) < Q{Ul)- That is, Ur is OC, and Q{Ur) 
is not greater than Q(Ul). The Riemann problem is solved by a combination of a backward 
traveling shock wave and a standing wave, with an intermediate state C/i as shown in Figure 
[H In this case, q = Q{Ur). 



Figure 16: An example for wave solutions to Riemann problem of Type 3 

Type 4 When Ur is in Region 4, pR < T{€r), Q{Ur) < Q{Ul), and Q{eR,r{eR)) < Q{Ul)- That is, 
Ur is UC, and the capacity at cr is smaller than Q(Ul)- The Riemann problem is solved by 
a combination of a backward traveling shock wave, a standing wave, and a forward traveling 
rarefaction wave, with two intermediate states Ui and U 2 as shown in Figure [T71 where U 2 = 
(e/j, r(efl;)). In this case, q = Q{eR,T(eR)), which is the capacity flow associated with cr. 

Type 5 When Ur is in Region 5, pR > r(e_R) and Q{Ur) < Q{Ul). That is, Ur is OC, and Q{Ur) 
is not greater than Q{Ur). The Riemann problem is solved by a combination of a backward 
traveling shock wave and a standing wave, with an intermediate state 17i as shown in Figure 
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Figure 17: An example for wave solutions to Riemann problem of Type 4 


m In this case, q = Q(Ur). 

Type 6 When Uh is in Region 6, pR > r(e/{) and Q{Ul) < Q{Ur) < (5(eL,r(eL)). That is, Ur is 
OC, and Q{Ur) is between Q{Ul) and the capacity at cr. The Riemann problem is solved 
by a combination of a backward traveling rarefaction wave and a standing wave, with an 
intermediate state Ui as shown in Figure fT^ In this case, q = Q{Ur). 

Type 7 When Ur is in Region 7, pR < T{eR), Q{Ur) < QicR, r(eL)), but Q{eR, T{€r)) > Q{eL, T{€l)). 
That is, Ur is UC, and the capacity at cr is between Q{Ur) and the capacity at cr. The 
Riemann problem is solved by a combination of a backward traveling rarefaction wave, a 
standing wave, and a forward traveling rarefaction wave, with two intermediate state Ui and 
U 2 as shown in Figure EiH where Ui = (eL,F(ei)). In this case, q = Q{U 2 ) = Q{Ui) = 
Q(ei,F(ei)), the capacity flow when e = tR. 
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Figure 18: An example for wave solutions to Riemann problem of Type 5 

Type 8 When Ur is in Region 8, Q{Ur) > (5(eL, r(eL)). The Riemann problem is solved by a combina¬ 
tion of a backward traveling rarefaction wave, a standing wave, and a forward traveling shock 
wave, with two intermediate state Ui and U 2 as shown in Figure HU where Ui = (eL,r(ei)). 
In this case, q = Q{U 2 ) = Q{Ui) = Q(eL,r(eL)). 

Type 9 When Ur is in Region 9, pR < r(e_R), Q{Ur) < Q{Ul), and Q{eR,T{eR)) < Q{Ul). That is, 
Ur is UC, and both Q{Ur) and the capacity at er are smaller than Q{Ul)- The Riemann 
problem is solved by a combination of a backward traveling shock wave, a standing wave, and a 
forward traveling rarefaction wave, with two intermediate state Ui and U 2 as shown in Figure 
HU where U 2 = {eR,T{eR)). In this case, q = Q{Ui) = < 5 ( 1 / 2 ) = <5(efi,r(efl)). 

Type 10 When Ur is in Region 10, pR < r(efi), and Q{Ul) < Q{eR,T{eR)) < (5(eL,r(pi)). That is, 
Ur is UC, and the capacity at cr is between Q{Ur) and the capacity at e^. The Riemann 
problem is solved by a combination of a backward traveling rarefaction wave, a standing wave. 
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Figure 19: An example for wave solutions to Riemann problem of Type 6 

and a forward traveling rarefaction wave, with two intermediate state Ui and U 2 as shown in 
Figure [23l where U 2 = (efl,r(e_R)). In this case, q = Q{Ui) = Q{U 2 ) = Q(e_R, F(eH))- 
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Figure 20: An example for wave solutions to Riemann problem of Type 7 



Figure 21: An example for wave solutions to Riemann problem of Type 8 
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Figure 22: An example for wave solutions to Riemann problem of Type 9 



Figure 23: An example for wave solutions to Riemann problem of Type 10 
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